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Abstract 

In this paper we show that the conventional diffusive shock acceleration mechanism 
for cosmic rays associated with relativistic astrophysical shocks in active galactic nuclei 
(AGNs) has severe difficulties to explain the highest energy cosmic ray events. We show 
that protons above around 2 x 10 20 eV could have marginally been produced by this 
mechanism in an AGN or a rich galaxy cluster not further away than around 100 Mpc. 
However, for the highest energy Fly's Eye and Yakutsk events this is inconsistent 
with the observed arrival directions. Galactic and intergalactic magnetic fields appear 
unable to alter the direction of such energetic particles by more than a few degrees. We 
also discuss some other options for these events associated with relativistic particles 
including pulsar acceleration of high Z nuclei. At the present stage of knowledge the 
concept of topological defects left over from the early universe as the source for such 
events appears to be a promising option. Such sources are discussed and possible tests 
of this hypothesis are proposed. 
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1 Introduction 



It is generally believed that cosmic rays with energies up to the "ankle" at around 3 x 10 18 eV 
are predominantly of galactic origin [fl] and that energies up to around 10 14 eV can be 
achieved by first order Fermi acceleration in shocks produced by supernovae exploding into 
the interstellar medium || . Recently the Fly's Eye detector revealed a change in the cosmic 
ray composition which is correlated with a dip in the total energy spectrum |1[ located at the 
ankle. Around this dip the spectrum first steepens and then flattens again to a spectral index 
of around 2.7 which is even smaller than the index of 3 corresponding to the spectrum below 
the steepening. The data are consistent with a superposition of a steep power law spectrum 
of heavy nuclei and a flatter spectrum of protons which overtakes the former component at 
energies above the ankle. It is expected that this latter high energy proton component is of 
extragalactic origin. Furthermore, on 15 October 1991, the Fly's Eye observed an event at 
(3.2 ±0.6) x 10 20 eV (la errors) 0, which is the event with the highest energy ever recorded. 
Interesting enough, the world's second highest energy air shower of (1.1 ±0.4) x 10 20 eV was 
recorded at Yakutsk || [7j located within 7.8 degrees from the Fly's Eye event (see Fig. 2). 
In this paper we will assume that these events were caused by relativistic particles. 

There is plenty observational evidence that AGNs and radiogalaxies contain relativistic 
termination shocks which are likely to produce high energy cosmic rays. It therefore seems 
natural to extend the standard diffusive shock acceleration scenario which works well for 
supernova shocks at lower cosmic ray energies to larger extragalactic shock scales in order 
to explain the origin of this higher energetic extragalactic component. However, it turns out 
to be hard to explain the highest energy events by this mechanism. An interesting option is 
provided by decaying or annihilating topological defects which could be left over from phase 
transitions in the early universe at temperatures corresponding to some Grand Unification 
scale. The rest of the paper is organized as follows: We first reconsider in section 2 the 
source spectrum cutoff energy E c for shock accelerated cosmic rays as a function of shock 
size and average magnetic field strenght on this lenght scale and evaluate it for some typical 
observational numbers. In section 3 we discuss propagation effects on protons and show that 
at least the 3.2 x 10 20 eV event is difficult to reconcile with the observational knowledge of 
typical extragalactic shock parameters in this acceleration scenario. We therefore discuss in 
section 4 some other options for these highest energy events. In section 5 we suggest that 
such events could alternatively be produced by topological defects. Finally we summarize 
our findings in section 6. 

2 The Source Spectrum Cutoff for Extragalactic 
Shock Acceleration 

In relativistic shocks the cutoff energy E c for the source spectrum of accelerated cosmic rays 
is in the test particle approximation always given by ZeBR, the product of the charge Ze 
of the cosmic ray particle, the magnetic field B and the size R of the shock, multiplied by 
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some factor of order unity || [UJ (we use natural units, i.e. c = h/2n = 1, throughout this 
paper). However, it turns out that for the highest energies the mean free path of the particle 
becomes comparable to the shock size R, which sometimes is not properly accounted for. 
We therefore calculate here our own approximation for E c . 



2.1 The Source Energy Cutoff 

The acceleration of a cosmic ray particle of energy E in an astrophysical shock is governed 
by the equation 

- = — , (1) 

dt T acc 

where T acc is the energy dependent acceleration time. In the statistical point of view the 
slope q of the energy spectrum dN/ dE oc E~ q of the particle flux is related to T acc and T esc , 
the mean (in general also energy dependent) escape time by 0. |9] 

9 = 1 + ^- (2) 



■ esc 



For first order Fermi acceleration at nonrelativistic shocks caused by supernovae, T acc is 
usually given by 

3 (D 1 D 2 



Tacc = — + — , (3) 

U\ — U 2 V U\ U 2 J 

where Ui, u 2 are the up- and downstream velocities of the shock and Di and D 2 are the 
corresponding diffusion coefficients, respectively. Diffusion is dominated by magnetic pitch 



angle scattering caused by inhomogeneities in the magnetic field [11 , [12[ . Therefore, the mean 
free path A is bounded from below by some multiple g of the gyroradius tl = E/(ZeB) and 
D\ and D 2 can for ultrarelativistic particles be estimated by 

For nonrelativistic shocks, g is usually set equal to 1 However, as we deal with 

the highest energetic extragalactic cosmic ray component, we have to consider relativistic 
shocks because they provide the most powerful accelerators. Monte-Carlo simulations of 
such relativistic shocks yield g ~ 40 |[D|, |13"| . Furthermore, the acceleration turns out to be 
enhanced compared to eq. @ by about a factor of 10 in highly inclined [[L3[] and by about a 
factor of 13.5 in parallel [jrjj relativistic shocks, respectively. 

Putting everything together and maximizing T acc from eq. ([3D we arrive at 

T > 9 E 
acc ~ 2.25 ZeB ' { ) 

On the other hand, as long as the diffusion approximation is valid, i.e. as long as A < R 
corresponding to E < Ediff = ZeBR/b, the escape time is given by T esc = R 2 / A, whereas for 
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E > Ediff the particles are freely streaming out of the shock region to a good approximation 
and T esc = R. Using eqs. (@) and @, we thus get 

Defining the cutoff energy E c as the energy where the source spectral index becomes 3 
(remember that the slope of the energy spectrum observed at the earth was around 2.7 in 
the region of highest energies) yields 



£ c = £ g=3 ~10 17 eVZ rM -77 • (7) 




This is compatible or even higher as compared to similar estimates 0, [TU], [T5], |T^, |T7|. We 
have assumed here that the magnetic field is parallel to the shock normal. If that is not the 
case there will be an electric field E = u x B in the shock rest frame (u is the shock velocity 
in the lab frame). This causes drift acceleration of charged particles to a maximal energy 
given by 



E max = ZeuBR ~ 10 18 eV Zu 



R \ B 



kpc J \/iG J 

which is around one order of magnitude larger than eq. (^) if u approaches the speed of light. 
However, the electric field E is expected to be much smaller in general due to plasma effects 
so that rather special conditions have to be fulfilled in order that such high energies can be 
approached. 

Throughout the rest of this section and the next section we will restrict our discussion to 
protons (Z = 1). We will comment on nuclei as possible candidates for events with energies 
above 10 20 eV in section 4.3. 



2.2 The Cutoff in Numbers 

Let us now look at some observational numbers for R and B and evaluate the corresponding 
cutoff energy for a proton. Cesarsky || cited the example of a galaxy encounter (NGC 
4038/39 |[L8|| ) as a location of a strong relativistic shock with a magnetic field of about 40 /iG 
on a scale of about 2 kpc, leading to E c ~ 8 x 10 18 eV, significantly too small to explain the 
origin of cosmic rays with energies as high as the Fly's Eye event of 3 x 10 20 eV. Potentially 
more interesting candidates for acceleration beyond 10 20 eV are revealed by the "hot spots" 
and radio lobes of CygA with B ~ 400 /iG, R ~ 3 kpc and B ~ 4 /iG, R ~ 300 kpc, 
respectively cited by Quenby []TU| which lead to a cutoff energy E c < 5 x 10 20 eV. There 



are also some indirect indications from gamma ray astronomy that in some quasars protons 
could be shock accelerated to energies of about 2 x 10 20 eV [GO]. 

There is one nearby object where pure application of formula (|7|) on observational indi- 
cations for R and B leads to an E c larger than 3 x 10 20 eV. This is if one takes the whole 
Virgo cluster with an extension of R ~ 10 Mpc and a intracluster magnetic field B ~ 1.5 /iG 
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which is compatible with observations |9|, [H]] and leads to E c ~ 1.5 x 10 21 eV. However, it 
is highly improbable that the whole Virgo cluster moving through intercluster space forms 
a relativistic shock of such an enormous extension, effective in coherent cosmic ray particle 
acceleration. 

Based on examples of the sort presented above there is some common belief that for 
protons the highest source energy achievable by diffusive shock acceleration in quasars and 
radiogalaxies is around 10 21 eV [|16 , 



J2|. We will nevertheless show in the next section 
that even if this is true, it is still difficult to explain the observed Fly's Eye event because of 
the information we have on its arrival direction. 



3 Propagation Effects on Protons above 100 EeV 
3.1 Overview 

Up to now we were only talking about the source energy spectrum. However, a proton trav- 
eling through space is in general subject to interactions, mainly with photons and magnetic 
fields. The latter effect leads to a curved path with a radius given by the gyroradius 




n = l-lpc TrnTzrf -77 ■ (9) 



The former effect leads to scattering and an effective energy loss as long as the proton energy 
lies above a kinematical threshold energy 

2e(l + cos 9) 

depending on the angle 9 between the incoming photon momentum and the negative nucleon 
momentum. Here is the nucleon mass, e is the typical energy of an incoming photon and 
Wth is the center of mass energy threshold for the particular reaction under consideration. 
The most important ones are electron pair production iV + 7— > N + e + + e~ with W^ e = 
niN + 2m e , and pion production, N + 7 — > N + n, with W£ h = rriN + m-*, with m e and 
the electron and pion mass, respectively. For a fraction of its propagation time depending 
on the effective neutron lifetime the proton will actually transform into a neutron which 
does, however, not have much influence on this treatment as protons and neutrons have 
similar interactions with the MBR. Thus, nucleons interact with photons of the microwave 
background radiation (MBR) producing e + e~ pairs and pions above an energy of about 
5 x 10 17 eV and 1.1 x 10 20 eV, respectively. The latter effect is known as the Greisen- 
Zatsepin-Kuz'min (GZK) effect [p3| |. There are also plenty of infrared and optical photons 



around luminous AGNs and galaxy clusters and especially in their central regions leading 
to correspondingly smaller nucleon threshold energies. However, their number density is 
in general not much bigger than that of the MBR photons and the cross section at the 
correspondingly higher center of mass energy is even smaller. This implies that possible 
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interactions with these higher energy photons does not substantially reduce the mean free 
path below 6 Mpc, the typical mean free path in the MBR f24}] . As this is much larger than 
typical galaxy sizes, we will restrict our considerations to interactions with the MBR. That 
means that we are on the safe side and get lower limits on energy losses and upper limits on 
the corresponding possible travel distances [ f25|j . 



3.2 Energy Losses 

We now ask the question in what distance range a source causing nucleon induced events of 
energy Eq on earth could be if the maximal source energy is E s . To this end we introduce the 
"longitudinal nucleon energy" E N = (1 + cos9)E N /2 > E th besides the nucleon energy E N 
in the comoving frame. The last inequality holds because of eq. ( |T0| ) as long as reactions are 
kinematically allowed. By performing a Lorentz transformation from the comoving frame 
to the center of mass frame corresponding to a gamma factor 7 cm ~ E^/W = En/{jti 2 n + 
46-Ejv) 1 / 2 one can see that after scattering the nucleon energy in the comoving frame is given 
by 

E 

E' N ~ lamiKm + Pern COS O ~ y^( E 'cm + Pern 9 cm ) , (11) 

where p' cm and E' cm are momentum and energy of the nucleon after scattering and 9 cm is 
the scattering angle, evaluated in the center of mass frame, respectively. Above 10 20 eV the 
energy loss is dominated by pion production []26] for which these quantities are related to 



W by (m 2 N + p'cm) 1 / 2 + (m 2 + p'cm) 1 / 2 = W. As we are only interested in an estimate we 
set En = -Eat/2 and neglect the energy dependence of the mean free path (which in the 
energy range we are interested in is a good approximation A ~ 6 Mpc in the following 
calculation. Because of eq. ( |TTD the energy change in a scattering event £(E, 9 cm ) relative to 
the lower energy is given by 

W,9 m ) = - + J cQs6 -1- (12) 

■^cm ' Pcm LAJO u cm 

Furthermore, we neglect energy loss due to cosmological redshift as we are dealing with non- 
cosmological distances at these high energies. By integrating from lower to higher energies 
one can show that the mean / and the variance A/ 2 of the distance as a function of E and 
E s can be estimated by 



cm , 



The first term in the variance is due to the fluctuation of the number of scatterings and 
the second one is due to the fluctuation of the energy transfer £ of eq. ( |T2| ) around its 
mean £(E) ~ W/E' cm — 1 averaged over the center of mass scattering angle 9 cm . We have 
numerically integrated eq. ([13]). The results are shown in Fig. la and lb. We see that for 
£o = 1.7 x 10 20 eV, the lowest possible energy for the Fly's Eye event the distance must be 
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smaller than ~ 100 Mpc and ~ 130 Mpc on the 3a level for a source energy E s = 10 21 eV 
and E s = 10 22 eV, respectively. For the best fit energy of 3.2 x 10 20 eV we get ~ 60 Mpc and 
~ 90 Mpc for the corresponding 3cr upper limits of distance for the same source energies. 

Thus from the energy point of view an AGN or a galaxy cluster constituting a large 
scale shock with the intercluster medium could be marginally able to cause events like the 
3 x 10 20 eV Fly's Eye event by shock accelerated protons, if it is not much further away than 
100 Mpc. Note that this number actually means the path length for which the distance is a 
lower limit which can be overtaken if the path is curved. However, we now show that the 
arrival direction of such a nucleon would then have to lie within about 10 degrees of the 
direction of its source if conventional wisdom about magnetic fields is used. 



3.3 Deviations from Rectilinear Propagation 

Let us first discuss deflection caused by magnetic fields. Unfortunately, not much is known 
about extragalactic magnetic fields. Faraday rotation measurements of extragalactic radio 
sources seem to suggest fields of the order of 10~ 9 G which could be homogeneous on large 
scales p7 |. Most estimates are of this order or below | 2~8|| , a more recent one being as low as 3x 



10~ n G [p9|1 . The bending angle a in radian for a proton traveling in a magnetic field satisfies 
a < / dl/ri where dl is the differential path length and tl is given by eq. (Qj. A 10~ 9 G 
field leads to a maximal bending angle of around 10 degrees for a proton with arrival energy 
of 3 x 10 20 eV. This maximum can only be reached if the magnetic field is perpendicular to 
the proton path and does not change its polarization considerably on a scale of ~ 100 Mpc. 
Otherwise the bending angle is reduced at least by a factor (dc/lOOMpc) 1 / 2 where d c is the 
magnetic field coherence length scale. Therefore, even if the typical intercluster field would 
be as high as 10 -8 G |30] the bending angle would still not be larger than 10° if the coherence 



length scale is ~ 1 Mpc. 

A proton can also be deflected by our own galactic magnetic field which is of the order 
of 3 fiG [E7, 31]. Based on radio telescope observations of Faraday rotations [31] its coherent 



component is supposed to have a cylindrical structure of diameter ~ 25kpc and height of 
order 1 kpc and being polarized in the direction of decreasing galactic longitude in the outer 
region. The random component is supposed to be of the same order of magnitude. Given 
the arrival direction of the Fly's Eye and the Yakutsk events shown in Fig. 2 we see that 
the corresponding path length through the field is less than 3 kpc resulting in a maximum 
bending angle of ~ 2°. 

We finally remark that even though magnetic fields in galaxy clusters as high as 10~ 6 G 
as already mentioned ]^T] could cause a significant deflection this does not influence our 
argumentation. It only means that if our proton encountered such "magnetic lenses" it 
should point approximately to the last encountered. There still has to be a nearby galaxy 
cluster in its arrival direction. Furthermore, the above mentioned fact that a proton partly 
transforms into a neutron during propagation even tends to decrease the bending further. 

The propagation direction could in principle also be changed by scattering with photons, 
but the scattering angles involved are much too small. To see that we use our notation from 
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above and note that the scattering angle 6 SC in the comoving frame obeys 

Maximizing with respect to 9 cm and using p' cm < (W 2 — m^) 1 / 2 < (AeE^) 1 ! 2 leads to an 
estimate independent from the final state X in the reaction iV + 7 — > iV + X, 



tan 6 sr . ^ 




iifi + ^l <io-, (15) 



where we have used e ~ 3 x 10 4 eV and E N > 3 x 10 20 eV. Because the mean free path 
at these energies is of order 6 Mpc |24| we expect only a few scattering events during trav- 



eling over a distance of 100 Mpc which thus never can lead to a significant change in the 
propagation direction. 

The arrival direction of the Fly's Eye event is given by a = 85.2° ± 1° and 5 = 48° ± 10°. 
No potentially interesting object with a sufficiently powerful shock acceleration engine of 
the scale discussed in the previous section is located within 100 Mpc in that direction. This 
can be seen from Fig. 2 where we show the directions to important nearby galactic objects, 
galaxy clusters and AGNs. One of the three prominent FR-II radio galaxies listed in [|T7|| 
and thought to contribute to the proton spectrum between 10 18 eV and 10 20 eV (i.e. below 
the GZK threshold), namely 3C111, has the coordinates a = 63.75°, 5 = 37.9° and lies 
thus (18.6 ± 5.1)° away from that direction. But it is at least 140 Mpc away. The only 
quasars within around 10° of the arrival direction are 3C147 and 3C159 P2"| , both of them 
at a distance of at least 900 Mpc. The only Seyfert galaxy within 10° is MCG 08.11 at 
an angular distance of (2.35 ± 8.4)° and a distance of at least 60 Mpc ||52" |. Although this 



is much nearer than the before mentioned quasars it produces a much lower radio flux at 
the earth and therefore seems also not likely to produce high cosmic ray fluxes. Under the 
quite improbable assumptions discussed at the end of the previous section the Virgo cluster 
seems to be the nearest possible candidate. However, it is located around 90° away from 
that direction. 



4 Other Options for Ultrahigh Energy Events 

We saw in the previous section that especially the Fly's Eye event is difficult to explain as a 
proton within the standard shock acceleration scenario. We therefore now like to discuss some 
other options beginning with secondary particles produced by shock accelerated protons. 

4.1 Secondary Photons 

Is the shower development of the highest energy Fly's Eye event consistent with what would 
be expected if it was caused by a photon? At these high energies the photon begins to 
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interact with the earth's magnetic field already above the atmosphere. Fitting the shower 



shape with a three parameter Gaisser-Hillas shower development function |33| gives a depth 
of first interaction of —100 gem -2 and thus seems to indicate a first interaction above the 
atmosphere. However, it could also indicate that the fitting function is simply inappropriate 
at these high energies. Indeed, it is now believed that fitting heavy nuclei induced showers 
can lead to similar negative values ||. Taking the LPM effect into account the average 
shower maximum is expected to be somewhat larger than for a proton induced shower with 
an average width about twice as large as the corresponding proton profile width [E3]. As 



the fluctuations are expected to be large this could still be compatible with the reported 
maximum at (815 ± 50) gem -2 . 

Photons of such a high energy have a secondary origin in the standard scenario as they 
are produced by decay of pions or e + e~ interactions which in turn are produced by the 
interactions of the cosmic ray protons with the MBR. The photon mean free path becomes 
comparable with or larger than that of the protons above a few 10 19 eV [ |35[| . The exact 
value of the photon to proton ratio depends on the universal radio background and the 
intergalactic magnetic field strength. The former leads to additional losses due to electron 
pair production. The latter leads to an inefficient electromagnetic cascade development 
7 + 7b — > e + e~ , e + e~ + jf, — > e + e~ + 7 (7^ is the background photon) due to synchrotron 
cooling of the electrons even for fields as low as 10~ 10 G. Typical estimates for the photon to 
proton ratios are considerably smaller than 1 above 10 20 eV |3(| |37]] . That makes the problem 
even harder as there have to be (even more abundant and of higher energy) primary protons 



acting as the source for such photons. Wolfendale p2| claims that the photon to proton ratio 
could exceed unity above ~ 3 x 10 20 eV if the source energy cutoff E c is much above 10 21 eV. 
This possibility thus runs into trouble with our discussion of E c in section 2. 

4.2 Secondary Neutrinos 

Could the highest energy Fly's Eye event have been an extragalactic neutrino produced as a 
secondary of a shock accelerated proton? Because neutrinos essentially lose no energy apart 
from redshift in going over cosmological distances it could have been produced by a proton 
interacting near its acceleration site thus avoiding excessive subsequent energy losses due to 
downscattering. However, it turns out that the neutrino yield above 10 18 eV is considerably 
smaller than one for all reasonable injection spectra. Furthermore, at the highest energies 
the spectral index observed at the earth is predicted to be 0.5 larger than the corresponding 



proton spectral index |j38[ . As even at these energies the neutrino nucleon cross section is still 
by a factor of at least 10 6 smaller than the nucleon nucleon (and also the gamma nucleon) 
cross section the average interaction depth is much larger than the atmospheric depth. All 
that leads to the conclusion that the event rate due to neutrinos should be much smaller 
than that due to protons at the same energy. One would also expect the neutrino induced 
showers to start predominantly at high depths i.e. near the horizon [j36fl . Therefore, if a 
neutrino caused the event then it was a very atypical one and we would expect much more 
proton and even photon events at the same energy. 
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4.3 Secondary Neutrons and Heavy Nuclei 

In the shock acceleration scenario neutrons are also produced as secondaries of protons or 
heavy nuclei as they are neutral and can not be accelerated directly. Furthermore, due to 
instability they have only a finite range which for 3 x 10 20 eV is about 3 Mpc. The puzzle is 
therefore not solved if our events are caused by neutrons. 

Looking at eqs. (|7|) and (0) one realizes that heavy nuclei can reach maximal energies 
which are higher by a factor Z compared to the protons. However, heavy nuclei lose energy 
not only due to the processes which dominate the energy loss of protons but also due to the 
giant dipole resonance which leads to photodisintegration. Above 10 20 eV the corresponding 
energy loss rates are about a factor 10 higher than those for protons [27)| and are typically 



due to proton stripping reactions. For example, a 56 Fe nucleus being launched with an energy 
of 10 21 eV will be below 10 20 eV after traveling 8 Mpc p4|. Thus, for our purposes nuclei are 



only interesting when they are of galactic origin. There are two galactic sites which could 
provide acceleration to interesting energies for heavy nuclei. The first one is the termination 
shock of the galactic wind caused by the milky way [^9|. This leads to maximal energies 
of ~ 3 x 10 17 eV 2(^600 kms" 1 ) 2 ^/!^ x 10 10 y)(-B/0.1 /j,G) § where u w is the galactic 
wind speed and Ta is the shock lifetime. Even for Z ~ 100 this is significantly too low 
and one is forced to use quite extreme parameters to reach beyond 10 20 eV. The second site 
would be even more natural to produce predominantly high energy heavy nuclei, namely 
young supernova remnants which form a pulsar wind shock ||. This is because pulsars can 
have quite large surface magnetic fields of order 10 12 G leading to fields of order 10 G on 
scales of 10 14 cm ||. The pulsar wind can be relativistic so that application of eq. (§) leads 
to E max ~ 10 18 eV Z which for Z ~ 50 is about one order of magnitude smaller than the 
highest observed energies. Adopting the more conservative estimate eq. (0) leads to a short 
fall of about two orders of magnitude. 

Indeed, as can be seen from Fig. 2, the Crab nebula lies near the arrival direction of our 
events and is expected to have magnetic fields of the order of 10~ 3 G on a scale of 2 pc ||0 
leading to estimates for the energy cutoff quite similar to the above mentioned. However, a 
calculation including the uncertainty of the latter one gives a relative angle of (25.9 ± 7.5)° 
which is more than 3a away. Furthermore, taking into account deflection effects due to the 
coherent galactic magnetic field component mentioned in section 3.3 increases the angular 
distance as the path should be bent towards the galactic north pole (see Fig. 2). For example, 
for Z = 5 the angular difference from the arrival direction would be (29.5 ± 6.9)°. It thus 
seems that only large bending by almost 360° could explain a possible origin from the Crab. 
For Z > 50 the Larmor radius is < 2kpc at these energies which indeed comes near the 
required amount of bending. In that case, however, the arrival direction is not expected to 
be correlated with the source location in a simple way and the propagation of these particles 
should better be considered as diffusion in the magnetic field. The source could therefore 
be any galactic site being able to produce the required source energy and this possibility 
can not be completely excluded although, as mentioned above, current models fall short in 
energy. It should be noted that the shower maximum at (815 ± 50) gem" 2 allows no definite 
distinction between a proton and a heavy nucleus induced shower as the expected numbers 
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for these options are ~ 850 gem -2 and ~ 775 gem -2 for iron, respectively |2], |3|. 

There are still problems left in interpreting the 3 x 10 20 eV event as caused by a heavy 
nucleus. As so it is possible that heavy nuclei could be disintegrated already at the source 
of acceleration p5| |. Furthermore, the Fly's Eye data between 10 18 eV and 10 20 eV suggest 



the transition to a lighter component as we already mentioned. 
4.4 Some Other Options 

There were some other suggestions how one could get to higher source energies. For example, 
Colgate ElJ claimed that in the relativistic plasma of AGN jets energies as high as 10 24 eV 



could be reached due to a plasma pinch effect similar than that used in Tokamaks. Note 
however that due to Fig. la and lb the source of the highest energy events could still not be 
much further away than 150 Mpc. Because the Larmor radius grows with energy the possible 
bending angle caused by magnetic fields could also not be enhanced significantly beyond 15 
degrees so the problem remains. 



It was suggested that high energy events could be caused by relativistic dust grains. 
The lateral shower profile caused by a dust grain entering the atmosphere is expected not to 
show a broad maximum but instead to have a more or less constant lateral spread as long 
as the grain remains large and energetic enough to produce secondary particles. However, 
the highest energy Fly's Eye event showed a quite "normal" shower development typical for 
a primary proton or possibly a photon which leaves the dust grain hypothesis to seem not 
very likely. Furthermore there is a tendency for these grains to break up by interactions with 
photons and gas atoms in the interstellar medium 



5 Cosmic Rays Produced by Topological Defects 

Topological defects (TDs) [[44j could have been formed in the early universe during phase 
transitions associated with spontaneous breaking of symmetries implemented in unified 
models of high energy interactions. Such TDs are magnetic monopoles, cosmic strings, 
domain walls, superconducting cosmic strings, textures, etc. TDs are topologically sta- 
ble but can nevertheless be destroyed due to physical processes like collapse or annihila- 
tion 1 45, 46, 47, . In that case the energy stored in the defects is released in the form 



of massive quanta of the fields like gauge fields and Higgs fields associated with the broken 
symmetry. These "X" particles released from the TDs would typically decay into quarks 
and leptons. Hadronization of the quarks would produce jets of hadrons containing mainly 
light mesons (pions) together with a small fraction (~ 3%) of nucleons. The gamma rays 
and neutrinos from the decay of the pions would thus be the dominant particles in the final 
decay products of the X particles. The mass mx of the X particles is typically of the order 
of the symmetry-breaking scale which in Grand Unified Theories (GUTs) can be ~ 10 25 eV, 
or even the Planck scale ~ 10 28 eV. The decay of the X particles released from TDs can thus 
give rise to nucleons, gamma rays and neutrinos with energies up to ~ mx, very much higher 
than what can be achieved by astrophysical shock acceleration mechanism. The cosmic ray 
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particles can thus be produced directly in this scenario, and no acceleration mechanism is 
needed. 

The production spectra of the nucleons, gamma rays and neutrinos in the TD scenario 
are determined by the physics of fragmentation of quarks into hadrons. Extrapolation j|5 



of QCD based hadronization models (which describe well the GeV scale collider data) to 
the extremely high energies gives a power-law approximation [47 , |48| , f4"9 [ to the differential 
production spectra with a power-law index q ~1.32 for nucleons as well as pions. The 
decay of the neutral pions thus gives a differential gamma ray production spectrum also 
with 5=1.32. It is to be emphasized, however, that there is a great deal of uncertainty in 
extrapolating the low energy QCD models of hadronization to the extremely high energies 
involved in the present situation. Moreover, the gamma ray production spectrum can be 
somewhat different from the proton production spectrum if one considers the gamma rays 
generated by the charged leptons (electrons and positrons) in the primary decay products 
of the X particles. The electrons and positrons coming from the decay of the charged pions 
in the hadronic jets also contribute to the overall primary gamma rays. The main point, 
however, is that the production spectra of cosmic ray particles in the TD scenario can in 
principle be considerably flatter than in the standard shock acceleration scenario. The latter, 
to recall, by and large produces differential production spectra with q > 2. 

One consequence of a relatively flat production spectrum in the TD scenario would be 
the "recovery" [49j of the evolved proton spectrum after the GZK "cutoff" ||23|[ . While this is 



heartening from the point of view of prospects for detecting protons above the GZK "cutoff" , 
too flat a proton spectrum may cause problems in that it may give rise to excessive gamma 
ray flux at much lower energies, as discussed below. In any case, as first discussed in Ref. [ 071 , 
the photon-to-proton ratio in the evolved spectra can be considerably larger than 1 above 
10 20 eV in the TD scenario J57J (because of the primary gamma rays which outnumber the 
protons by a factor of at least 10 at production, and also because of higher transparency of 
gamma rays relative to the protons at these energies), and so the cosmic rays above 10 20 eV 
are predicted to be mainly primary gamma rays rather than protons. 

Gamma rays as well as protons of ultrahigh energies generate lower energy gamma rays 
by 7 — lb and p — % collisions with the photons (7^) of the background radiation fields. 
The electromagnetic component of the energy lost by the photons and protons in these 
collisions cascades down to lower energies by electromagnetic cascading in the universal radio 
background (URB), the microwave background (CMBR), and in the infrared background 
(IRB) (in order of decreasing energy of the propagating photon). Recently it has been 



realised pUL 51, B2[ that the measured flux of extragalactic gamma rays in the 100 Me V 



region [53|| provides an upper limit on the total energy density of the cascade-initiating 



electromagnetic radiation that can possibly be released in the universe due to p — 7^ and 
7 — 7b interactions. This, in turn, restricts the shapes of the proton as well as the primary 
photon spectra in the highest energy region. The authors of Ref. [^0] claim that a proton 
spectrum with q = 1.32 at injection and extending to 10 24 eV would by itself give rise to a 
100 MeV gamma ray flux exceeding the measured flux by a factor of 2 if the evolved proton 
spectrum is normalized [37, |49j to observed particle flux at 4 x 10 19 eV. This question has 
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recently been studied in detail by a careful numerical calculation of the cascading process 
including the gamma rays generated by both p — 7^ as well as the 7 — 76 processes. It is 
found that whether or not the predicted 100 MeV gamma ray flux exceeds the measured 
value depends strongly on the level of the IRB as well as on its cosmological evolutionary 
history both of which are rather uncertain, and so a firm conclusion in this regard cannot be 
drawn at this stage. Nevertheless, the authors of Ref. |51| have suggested that the possible 
problem arising from requirement of consistency with the measured 100 MeV gamma rays 
can be avoided if the cosmic rays above 10 19 eV are mainly gamma rays and not protons. The 
preliminary analysis of Ref. [fy] shows that this is possible provided the injection spectrum 
of the primary gamma rays above 10 20 eV in the TD model is made somewhat steeper 
(q ~ 2.4) compared to the protons (q ~ 1.5) and 7/p is demanded to be ~ 60 at injection 
so that the proton component is made negligible compared to the photons. The average 
multiplicity in the hadronic jets arising out of the decay of the X particles is also required 
to be somewhat higher than what naive extrapolation of the low energy QCD based models 
of jet fragmentation indicates. While all these phenomenological requirements need to be 
substantiated on more theoretical grounds, the general conclusion that seems to arise from 
the above discussion is that the highest energy cosmic ray particles in the TD scenario should 
be mainly gamma rays and not protons. And, of course, primary neutrinos [f|^J should be at 
least as abundant as the gamma rays, perhaps even more. 

Could the 3 x 10 20 eV Fly's Eye event be a primary gamma ray due to TD collapse or 
annihilation? As already mentioned in Section 4.1 above, the shower development is not in 
contradiction with what is expected for a primary gamma. The resulting electromagnetic 
shower can in fact be very similar [K3] to a proton-induced shower, although some differential 
parameters, e.g., muon/electron ratio at large distances from the core of the shower can in 
principle be used for effective separation [^] of these photon-induced showers from the 
proton-induced ones. 

How could one distinguish between the TD option and the galactic heavy nuclei hypoth- 
esis which seems to be the least problematic option within the standard picture? Heavy 
nuclei are expected to produce substantially more muons compared to gammas of the same 
energy 0. It should therefore also be possible to draw a decision between these options as 
more statistics is available at these highest energies. 

The lack of any obviously identifiable astrophysical source for the event is not a problem 
for the TD scenario because TDs are not necessarily expected to be associated with any 
astrophysical sources such as galaxy clusters or AGNs. The TD model thus seems to offer 
an attractive option in this regard. 

It is, however, expected that the same TD annihilation event would also produce lower 
energy gamma rays which would arrive at earth at roughly the same time and with same 
arrival direction as the 3 x 10 20 eV Fly's Eye event. Unfortunately, the CASA array ]55] 
capable of detecting such gamma rays was not operating at the time when the above Fly's 
Eye event was recorded. However, the CYGNUS array [5fJ capable of detecting gamma rays 
above about 100 TeV was operating and it detected no event []57| that can be associated 
with the Fly's Eye event. If the integral primary gamma ray spectrum between 10 14 eV and 
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10 20 eV due to TD annihilation is taken to be approximately proportional to E a , (a > 0), 
and if one (optimistically) takes the flux at 10 20 eV as ~ 1 per 1000 km 2 , then above 10 14 eV 
and in an area of 0.02 km (roughly the area of the CYGNUS array) one would expect an 
integral flux F(E > 10 14 eV) ~ 2 x 10 60-5 events per 0.02 km 2 . The non-detection by the 
CYGNUS array of any gamma ray event in the 100 TeV region coincident with the Fly's Eye 
event can then be interpreted in terms a being < 0.78 (i.e., a relatively flat spectrum) in 
the TD model at energies below 10 20 eV. (Note that in conventional scenarios a is usually 
taken to be ~ 1). For example, if in the TD model one takes a ~ 0.32 |4£j and neglects 
the attenuation due to interaction with the CMBR (thus overestimating the expected flux), 
then F(E > 10 14 eV) ~ 0.0017 events in 0.02 km 2 , and so CYGNUS may have missed the 
event. This point, however, needs further investigation, and will be discussed elsewhere. 



6 Conclusions 

We are lead to the conclusion that protons arriving at the earth with energies of 3 x 10 20 eV 
or above are very likely to have come from an AGN or a galaxy cluster not further away 
than 100 Mpc if they were produced there via the standard diffusive shock acceleration 
mechanism. Even then the necessary conditions to be fulfilled in such relativistic strong 
shocks seem highly improbable as long as the shock parameters have to be compatible with 
observational data. Furthermore, the arrival direction of such protons have to be within 
around 15 degrees in the direction of their source. The 3 x 10 20 eV Fly's Eye event and the 
highest energy Yakutsk event were therefore very likely not protons produced within the 
standard diffusive shock acceleration scenario as they do not point to some possible source 
being nearer than 100 Mpc. 

Some other explanations for such events like that being produced by secondaries of shock 
accelerated protons were discussed. Within the astrophysical shock scenario the most promis- 
ing, although also problematic option seem to be heavy nuclei of galactic origin which could 
be accelerated in pulsar wind shocks driven by young supernova remnants. 

We therefore conclude that at least some improvements in the understanding of the 
current acceleration picture have to be made in order to explain the highest energy cosmic 
rays observed. It seems possible that a completely new production mechanism for such 
particles is necessary. We suggested that the TD model could be a promising option. It is 
curious that such an exotic option seems to have less difficulties in explaining these ultra- 
high energy particles. All other current options appear to require suspension of belief in 
seemingly well substantiated observational numbers or indicate incomplete understanding of 
the underlying physical process. We list all the options discussed here together with their 
problems in Table 1. 
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protons and secondary 7's and z/'s 



Sources 



maximal energy 



shock acc. 



drift acc. 



additional 
problems 



AGNs 



distance+direction 



pulsars 



2 x 10 17 eV( 



RB 



0.002 pcG 



2 x 10 iS eV 



RB 



0.002 pcG 



direction 



galactic wind 



3 x 10 17 eV ( 600k " ns -i) (juJopiOy) (o.i^G 



galactic heavy nuclei 



Sources 



maximal energy 



shock acc. 



drift acc. 



additional 
problems 



pulsars 



2xl0 17 eVZ( 5 ^ 



2xl0 18 eVZ( o Jf^ 



direction 



galactic wind 



3 x 10 17 eV Z (eoo^g-i) (i. 5 >aoiOy) (o.i^g 



mainly 7's, some protons and z/'s 



maximal energy 



problem 



TD's 



> 10 24 eV 



exotic 



Table 1: Maximal energies and problems with the options discussed in this paper. 
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Figure Captions 



Figure 1A and IB: Distance of source versus source energy for protons of arrival energies 
of Eq = 1.7 x 10 20 eV and 3.2 x 10 20 eV, respectively. Plotted from bottom to top are the 
average distances and the maximal distances at the la, 2a and 3a level, respectively. 

Figure 2: Arrival directions of the highest energy events seen by Fly's Eye and the Yakutsk 
experiment in galactic coordinates [pD ]. This plot is centered around the galactic anticenter 



with the middle horizontal line being the projection of the galactic plane. Also shown are 
nearby galaxy clusters (big circles), AGN's (small circles) and galactic supernova remnants 
(light circles). As discussed in section 3.3 the galactic magnetic field is supposed to have 
a coherent component near the galactic plane which in the outer region is polarized in the 
direction of decreasing galactic longitude. Therefore, the apparent arrival direction of cosmic 
rays coming from one of the objects shown should be shifted by ~ 1.5° Z(3 x 10 20 eV/E) to 
lower latitude. 
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